N-dimensional alternate coined quantum walks from a dispersion relation perspective 
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Abstract 

We propose an alternative definition of an N-dimensional coined quantum walk by generalizing a recent proposal [Di Franco 
et at., Phys. Rev. Lett. 106, 080502 (2011)]. This N-dimensional alternate quantum walk, AQW^p^j, in contrast with the 
standard definition of the N-dimensional quantum walk, QW^j^j, requires only a coin-qubit. We discuss the quantum diffusion 
properties of AQW(2) and AQW(3) by analyzing their dispersion relations that reveal, in particular, the existence of diabolical 
points. This allows us to highlight interesting similarities with other well known physical phenomena. We also demonstrate 
that AQW(-3-) generates genuine multipartite entanglement. Finally we discuss the implementability of AQW^j^j. 

PACS numbers: 



In both its standard forms, the coined fl| and the con- 
tinuous one 2] , the quantum walk is the quantum version 
of a classical random process, described by the diffusion 
and the telegrapher's equations, respectively In the 
coined quantum walk - the process we consider here - 
there is a system (the walker) that undergoes a condi- 
tional displacement, to the right or left, depending on 
the output of a coin throw, as in the random walk. But 
differently from its classical counterpart, here both coin 
and walker are quantum in nature; hence the walker po- 
sition coherently spans over the lattice sites. The one- 
dimensional coined quantum walk - QW^j^^ for short ~ 
has been studied from many different perspectives, espe- 
cially from the quantum computational point of view [J] . 
In the last few years it has also received increasing ex- 
perimental attention 043 j includingthe case of quantum 
walks with more than one particle [8|. 

The situation is quite different when dealing with N- 
dimensional quantum walks, QW^^) for short. They were 
first discussed by Mackay et al, who introduced them in 
complete analogy with QW^^ (see also Ref. 
As defined in Ref. 0, QW^^^ requires the use of a 2N- 
dimensional qudit as coin, as well as a coin operator rep- 
resented by a 2N x 2N unitary matrix. This introduces 
increasing complexity in the process as N grows, espe- 
cially from the experimental viewpoint ll|, |l2| , but also 
from the theoretical one M\- However, Di Franco 
et al. [3] have recently proposed an alternative two- 
dimensional quantum walk, namely the alternate quan- 
tum walk - AQW for short - that is simpler than the 
standard one. In AQW the coin is a single qubit, as in 
QWj]^-), and each time step is divided into two halves: in 
the first one the coin is thrown (i.e., a Hadamard transfor- 
mation is applied on the coin-qubit) and the conditional 
displacement along x is consequently performed; then, in 
the second half of the time step, the coin is thrown again 
and the conditional displacement along y is performed. 
Hence, in AQW, the four-dimensional qudit of QW(2) is 



replaced by a single qubit, the price paid for that being 
to double the number of sub-steps per single time step. 
Quite unexpectedly, AQW reproduces the same spatial 
probability distributions of QW(2) when the Grover coin 
is used - Grover-QW(-2) for short - for a set of particular 
initial conditions, precisely those for which the charac- 
teristic localization of Grover-QW(-2) does not occur. In 
Refs. 1^ 3| analytical demonstrations of the (limited) 
equivalence between Grover-QW('2) and AQW are given. 

Here we generalize AQW to N-dimensions: We de- 
fine the alternate N-dimensional coined quantum walk 
- AQW(N) for short - as a quantum walk in which the 
time steps are divided into N sub-steps. In each of these 
sub-steps, the coin throw is followed by the conditional 
displacement along one of the N dimensions. In a single 
time step of AQW(ivi): the qubit-coin is therefore thrown 
N times, but this is clearly simpler than throwing a sin- 
gle 2N-dimensional coin, from both the experimental and 
analytical points of view, in particular for large N. We 
show below that, besides a simpler experimental imple- 
mentability as compared with QW^j^-j, AQW(]sr) has a 
rich dynamics that is also simpler to understand than 
that of QW(N-). We provide some analytical results con- 
cerning the evolution of the probability distribution for 
cases N = 2, 3, paying special attention to the processes' 
dispersion relations (that reveal the existence of diabol- 
ical points), and demonstrate the existence of genuine 
multipartite entanglement in AQW(3). We end this let- 
ter with a brief discussion concerning the experimental 
implementability of AQW(n)- 

In order to introduce AQW(n) formally, let 
represent the state of the system at (discrete) time t. 
Vector is defined in the compound Hilbert space 
Tic <8) 'Hp, where Tic and Hp are the Hilbert spaces 
for the coin-qubit and lattice sites, respectively. Being 
He and Hp spanned, respectively, hy {\c) ,c ~ u,d} and 
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{\x) , X e Z^}, we can write 



(1) 



\x;c) = \x) (E) |c), where cs^t is the probabihty amphtude 
for the walker to be at site x = (xi, . . . , xjv) at time t 
with the coin in state c. The probabihty of finding the 
walker at site x and time t is t = \us + Ids ■ 



The state evolves as |'0)t+i = ^{n) 



)j, with {/(AT) = 
D]\[CnDn-iCn-i ■ ■ ■ DiCi a unitary operator. Here, 
operator Ci is the coin operator^ acting only in He, whose 
more general form is 



Ci = cos I 



[\u) (t 



\d) {d\ 



(2) 



3in6', {e'"'\u} {d\-e'^^ \d} {u\) 



with (ai, Pi,9i) arbitrary reals and i — I, . . . , N (notice 
that, for ai — /3i — and Oi — 7r/4, Ci is just the 
Hadamard transformation). Operator Di is the condi- 
tional displacement operator along direction xi, which 
we write as 



D, = 



[\x + ni;u) (f ; -u| + |x - n^; d) (f ; d\] , (3) 



where fti is the unit vector along direction Xi 



Equation \^) 



t+i 



U, 



(TV) 



can be expressed as a two- 



dimensional map relating probability amplitudes cs,t+i 
with cs'^t with c, c = u,d and x' nearest neighbors to 
X. This map admits two plane-wave solutions of the 
form col = exp [i • x — cj'^^t)] where 

'fq,± = col{u^^±,d^^±) are time independent vectors, qis 
the pseudo-momentum, with qi € (— 7r,7r], and w*^^^ are 
two frequencies determined by the dispersion relation. 

The dispersion relation is most relevant because \'ip)^ is 
entirely determined by it, given \ip)f.^Q. Moreover, when 
the initial state extends over a finite set of points in the 
lattice, especially when it is modulated by a smooth func- 
tion of space, the dispersion relation is particularly useful 
for predicting the evolution of the initial wave-packet, 
due to the relatively well defined group velocity (given 
by the local gradient of the dispersion relation curve). 
In this case, long-wavelength continuous approximations 
are very well suited and useful for envisaging the long 
time behavior of the probability distribution. This has 
been discussed in detail for QWj^^ in Ref. [Si (see also 
Refs. i,li3l). Here, we will limit ourselves to a quali- 
tative discussion of what the dispersion relation suggests 
for cases N = 2, 3. 

Let us first consider AQW(2)- By proceeding as stated 
above, one obtains the following dispersion relation 

cos n — C1C2 cos (u + v) + S1S2 cos {u — v) , (4) 

where q — cosOi, Si — sin^^, u = qi + (/3i + /2, 
V = q2 + /32/2, and 17 = w - (^1 +^2 + ^2) /2. Notice 



that phase ai does not appear in Eq. ([4]), hence it is 
irrelevant. As for the phases that do appear, they just 
entail a translation of the frequency and spatial quasi- 
momentum. Only 9i and 62 are dynamically relevant 
parameters. 

In Fig. 1(a) we have represented the dispersion curves 
for 9i ~ 62 ~ 7r/4. The most relevant features are the 
existence of a number of saddle-points (for which the 
group velocity is zero) , together with regions of maximum 
slope, which equals 0.5 (hence the maximum velocity in 
AQW(2)) and, most importantly, the existence of diabol- 
ical points (DPs in the following). DPs are conical inter- 
sections involving a degeneracy [20[. Interestingly, when 
Si 7^02, the DPs disappear, as shown in Fig. 1(b). This 
suggests that the existence of DPs could be particularly 
sensitive to decoherence effects in the coin mechanism. 

The dynamics around the points of null and maximum 
slope can be analyzed as in Ref. [11] for QW^^), i.e., 
in terms of known solutions of simple linear wave equa- 
tions. One can then envisage more or less straightforward 
generalizations of the results there discussed to the two- 
dimensional case (we discuss this elsewhere for Grover- 
QW(2) [223 )• But the existence of DPs is particularly 
appealing and constitutes a qualitative difference with 
the one-dimensional case. This geometric object, the DP 
(that takes its name from the diabolo-like shape of the 
conical intersection), appears in physics in quite differ- 
ent contexts like, for instance, the electronic spectrum of 
polyatomic molecules , quantum triangular billiards 
[20|, conical refraction in crystal optics 21|, or the dis- 
persion relations for massless fermions (Dirac electrons) 

Hi], and the 
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in QED and for electrons in graphene 
diabolo is associated with some of the remarkable phe- 
nomena appearing in those systems. We can expect to 
find in the dynamics of AQW(2) some parallelism with 
phenomena present in the aforementioned systems. In 
this sense. Fig. 1(c) shows the propagation of an initially 
extended state (with a Gaussian probability distribution 
of uhwhm — 7) after 90 time steps for an initial coin 
state col(l,i) /\/2 equal for all the populated sites. A 
homogeneous ballistic propagation from the origin with 
perfect circular symmetry is clearly visible, which recalls 
the ring shape of the beam exiting a biaxial crystal un- 
der conical refraction conditions [21]. It is interesting to 
compare this dynamics with that shown in Fig. 1(d), 
which has been obtained for the same initial conditions 
but with 9i = 7r/4 ^ 82 = n/3. In this case the diabolo 
is lost and the branches of the dispersion relation show 
a parabolic shape [see Fig. 1(b)] that leads to evolutions 
typical of linear optical diffraction [11] . It is interesting 
to mention that a similar controlled disappearance of the 
DP has been experimentally observed in graphene 24]. 

Before moving to a higher dimension, let us revisit 
the relation between AQW(2) and Grover-QW(-2) from 
the dispersion relation perspective. The dispersion re- 
lations for Grover-QW(2) consist of four sheets, because 



2 



the coin space is four-dimensional, and can be found in 
Refs. IJ, |22|. In our notation, they read uji,2 = 0,7r and 
u!3^4 = ± arccos [(cos u + cos v) /2]. Remarkably, the two 
sheets uj^a coincide with those of AQW(2), Eq. (jH), for 
01 = O2 ^ 7r/4 after identifying {u,v) in Grover-QWj-j) 
with {u + v,u — v) in AQW(2), i-C, both dispersion re- 
lations coincide for these parameters up to a 7r/4 rota- 
tion of the pseudo-momentum. The two other roots in 
Grover-QW(2), wi,2 — 0,7r, are constant, which means 
that the projections of the initial state onto the corre- 
sponding eigenvectors will not evolve in time. This is the 
origin of localization in Grover-QWr^-i for most initial 
coin states, as already noticed in Ref. We conclude 
that, whenever the initial state in Grover-QW(-2) does not 
project onto the eigenvectors governed by uji,2, Grover- 
QW(2) and AQW(2) are isomorphous for di = 62 = 7r/4. 
This is our proof of the (partial) equivalence between the 
two versions of the process. 

Let us now move to AQW(3) . The dispersion relation 
is governed by 



sin Q, = ci [C2C3 sin (u + w + w) + S2S3 sin (u — w + w)] 
+ si [C2S3 sin (u + V — w) — S2C3 sin (u ~ v ~ w)] , 

(5) 

where Ci — cosOi, Si — sinff^, and (u,w,u',r2) — 
(91, 92, 93, w) -I- {5qi,5q2,5qs,5uj). Here 25qi = ai + 
^2, 25(72 = a2 + (33, 2Sq3 = 03 + /3i, and Sio = 
— {Sqi + 6q2 + 5q3) ■ From ([5]) two dispersion relations 
are obtained, namely w^^-* — Q. and w*-"^ = tt — w'-^^ 
As in AQW(2), some phases in d are absent in the dis- 
persion relation (hence they are irrelevant) and the effect 
of the rest of phases in Ci is just a displacement of the 
dispersion relation surfaces in the {gi, (72, 93, w} space. 
Eq. ^ is simpler for 9^ = 7r/4, i = 1, 2, 3 

sin u) = [sin v cos [u — w) -\- cos v sin {u + w)] , (6) 
v2 

in which case there are eight degeneracies (occurring 



when 



,(+) 



uj^^^i = 7r/2) at {udp,vdp,wdp) = 
{(a, a, a) , (—a, —a, b) , (—a, b, —a) , (6, —a, —a)}, with 
a,b = 7r/4, 37r/4 and a ^b. These are the 3D equivalents 
of the DPs discussed for AQW(2). As in AQW(2), the 
DPs disappear when phases 9i are different. We mention 
that several crystallographic structures have recently 
been proposed for obtaining 3D-DPs (Dirac-semimetal 
in 3D) [il. 

In Figs. 2(a,b) we represent two bidimensional pro- 
jections of the probability distribution corresponding to 
the propagation of a walker initially localized at the ori- 
gin. The width of the distribution grows linearly with 
time, along all the three spatial dimensions, as it hap- 
pens for lower dimensionality. In order to show the effect 
of the DP, Figs. 2(c,d) show the same projections as 
Figs. 2(a,b) when the initial condition is not localized 



but extended. Again we choose a Gaussian distribution 
with ohwhm = 7, whose pseudo-momentum is centered 
at one of the DPs. We observe a symmetric ballistic dy- 
namics in the {xi,X2) plane. Fig. 2(c), that resembles 
that of Fig. 1(c). Notice the existence of two concen- 
tric bright rings, as it occurs in conical refraction (20| . 
However, the probability is not equally symmetric in the 
{xi,X3) plane. Fig. 2(d), which reveals an intrinsic lack 
of symmetry in AQW(3): indeed it can be shown that 
it is not possible to find any initial coin state that leads 
to a symmetric propagation in all directions (the cho- 
sen initial coin state, (0, 1) for all the initially populated 
sites, leads to a symmetric distribution in the planes with 
constant 0:3, but not in those for which X3 varies, as the 
figures show). 

One could wonder whether AQW(3) is a process similar 
to Grover-QW(-3-) as it happens for N = 2. The answer 
is negative: we have derived and compared the disper- 
sion relations of both processes for N = 3 and they are 
different. Hence N = 2 is a singular case in this respect. 

Having a richer and more complex structure than the 
bipartite case 26[, multipartite entanglement has re- 



cently attracted a lot of interest in the scientific commu- 
nity. Glearly, a feasible system able to generate a proper 
amount of genuine multipartite entanglement could be 
a valuable benchmark for this rapidly growing research 
field. Hence a most relevant question is whether quantum 
walks with N > 2 do exhibit genuine multipartite entan- 
glement or not. We have investigated this for N = 3. For 
analyzing its generation in AQW(3) , we have first to trace 
out the state of the coin. We are then left with a density 
matrix in the composite Hilbert space T-Lxi ® T~Lx2 ® 
(each subspace corresponding to a direction of the walk). 
We evaluate the multipartite entanglement present in 
this composite system by means of the tripartite neg- 



ativity [27[. This is defined as the geometric average of 



the three negativities that are obtained by considering 
the three possible bipartitions of the total system, giving 



-23-^2-13^^3-12- 



(7) 



Here Ni^jk is the negativity of the composite system 
{i, J, k} corresponding to the bipartition in the subsystem 
{i} and the subsystem {j, k}. Each Hilbert direction- 
subspace has a dimension growing with the number of 
time steps, so we use the generalization of the negativity 
for higher-dimensional systems (so as to have < iV < 1) 
[iij. We have evaluated A^^^^ in AQW(3) , with the walker 
starting at the origin and initial coin state col(l,i) /V^, 
for a number of time steps t up to 10, obtaining the plot 
in Fig. 2(e). Even if the number of time steps considered 
here is not so large (due to the dimension of the total 
Hilbert space, the computational power required for the 
evaluation grows rapidly with t), it is easy to check that 
AQW(3) is able to properly generate genuine multipartite 
entanglement. 
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Let us finally discuss the implementability of AQW(n) • 
Realizing QW^^) is quite demanding because of the com- 
plexity of implementing coin-operators to transform the 
needed 2N-dimensional qudit. On the other hand, in 
AQW(N): (i), a single coin-qubit is required indepen- 
dently of N; (ii), two-dimensional transformations of the 
qubit are easy to implement ^^-d (iii), the sequen- 

tial application of operators DjCj, j — I, . . . , N, makes 
the implementation of AQW(n) similar to that of QW^j^-, 
provided that all N dimensions could be multiplexed into 
a single one (similarly to what Schreiber et al. 12 1 



have recently done in their pioneering implementation of 
QW(2))- Indeed one could even implement AQW(n) with 
a constant number of physical elements independently of 
N if one has sufficient control on the experimental device. 

In order to illustrate the idea, we generalize the ideal- 
ized device already discussed in Refs. [H, [2^, which is 
similar to that actually used in Ref. @ . Consider a long 
enough optical cavity containing two electro-optic mod- 
ulators (EOM) whose roles are (i), EOMl performs the 
coin-operator C (i.e., makes a unitary transformation of 
the light polarization state, which plays the role of the 
coin-qubit in this implementation of the walk) ; and (ii) , 
E0M2 performs the conditional displacement Z), which 
consists in up/down shifting the carrier frequency of the 
light pulse depending on its polarization. The light pulses 
entering the cavity are assumed to be much shorter than 
the cavity length, and the frequency shifts introduced 
by E0M2 must be large enough for avoiding any fre- 
quency overlapping between pulses. With such a device 
not only QW(i) can be implemented, see [EES, US], also 
AQW(N) could be implemented for different values of N 
by properly programing the operation of the EOMs, with- 
out the need for additional elements. For example, in 
order to perform AQW ^2): the first half of the time step 
is implemented during one cavity round-trip of the light 
pulse (during which EOMl and E0M2 implement DiCi)] 
then, in the subsequent cavity round-trip, the settings of 
both EOMs are changed in order to perform a different 
coin operator and a different frequency displacement im- 
plementing 1)2(^2. Importantly, the frequency displace- 
ments in Ci and C2 must be different enough in order 
to multiplex a large number of steps [fij . Only technical 
limitations seem to limit the extension of the procedure 
to higher N. In any case, with this oversimplified discus- 
sion, we are not claiming that the device just outlined 
is the most appropriate for implementing AQW(n). In- 
deed a suitable modification of the fiexible device used 
by Schreiber et al. 0, [13 would probably be a more 
promising option. With our discussion we just want to 
emphasize that a single and conceptually simple device 
could implement alternate quantum walks with tunable 
dimensionality. It is much interesting that such a device 
would be a source of genuine multipartite entanglement. 

In conclusion, we have introduced the N-dimensional 
alternate quantum walk and discussed some of its proper- 



ties through the analysis of the dispersion relations that 
reveal, in particular, the existence of diabolical points. 
We have also demonstrated that for N = 3 the process 
generates genuine multipartite entanglement. We have 
finally discussed the implementability of the process that 
could be possible with physical resources that do not nec- 
essarily grow with the dimensionality of the walk. 
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Figure captions 

Fig. 1. 3D views of the two branches of (HJ, for 
6»i = 6I2 = 7r/4 (a) and 6*1 = 7r/4 and 02 = tt/3 (b). 
Propagation in AQW(2) after 90 time steps of an initial 
state with a Gaussian probability distribution of width 
<^HWHM = 7 and initial coin state col(l/\/2, i/V2) equal 
for all populated sites, for 9i =62 = 7r/4 (c) and di = 7r/4 
and 92^tt/3 (d). 

Fig. 2. Propagation in AQW(-3-) after 90 time steps 
of (a,b) a localized initial state, and (c,d) a spatially ex- 
tended initial state (Gaussian probability distribution of 
CHWHM = 7) for 61 = O2 = IT / ^ and initial coin state 
col(0,l). In (e) the tripartite negativity is shown as a 
function of the number of steps. 
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